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ABSTRACT 

"T^he  problems  of  processing  of  satellite  radiophysical 
measurement  and  mapping  of  ocean  surface  and  atmospheric  characte¬ 
ristics  in  a  monitoring  regime  are  discussed.  The  principles  of  two- 
dimensional  processing  of  multiparametric  and  multichannel  measure¬ 
ments  with  regard  to  the  requirements  of  operativeness  and  reliabi¬ 
lity  of  the  monitoring  system,  accuracy  of  solution  of  identifica¬ 
tion  and  detection  problems  are  considered.  The  structural  solutions 
on  organisation  of  a  real-time  data  processing  are  suggested. 

The  Importance  of  the  development  of  operative  methods  of  geophy- 
sioal  parameter  determination  in  a  monitoring  system  technology  is 
noted.  As  an  example  of  microwave  radiometric  data  processing,  we 
describe  the  computer  technique  used  on  the  satellite  Cosmos-1151 
for  measuring  the  atmospheric  and  oceanographic  parameters. 

The  estimations  of  accuracy  of  geophysical  parameter  determination 
are  given. v 

The  study  of  environment  in  global  scales  requires  the  organization 
of  global  monitoring  systems.  Such  systems  may  have  a  complex  hierarchy 
structure  Involving  a  satellite,  earth-based  receiving  and  data  processing 
complexes,  aeroplanes,  ships  etc.  The  operative  detection  and  identification 
of  anomalies  (pollutions,  fires,  typhoons,  hurricanes,  tropical  cyclones  on 
the  Earth's  surface  etc.)  are  regard  to  be  one  of  the  principle  purposes  of 
such  systems. 

The  problem  of  ditecting  anomalies  on  the  Earth's  surface  requires  for 
its  solution  the  application  of  up-to-date  automatic  decision  machines  for 

processing  the  large  arrays  of  satellite  measurements.  This  problem  involves 
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the  organization  of  in-line  data  prooeasing  at  all  detection  system  levels 
under  time  conditions  ensuring  an  uninterrupted  matching  of  operation  of  all 
cybernetic  devices  of  the  system  taking  into  account  their  functional  limita¬ 
tions  (computer  speed,  memory  capacity,  etc.) 

The  search  for  anomalies  on  the  Earth's  surface  involves,  on  the  whole, 
the  following  main  stages:  the  in-line  data  enters  the  computer  in  the  form 
of  m-dimensional  vector  T,  where  the  number  IT*  is  determined  by  the  number 
of  information  channels  and  by  the  number  of  measured  parameters.  Then  accor¬ 
ding  to  the  computer  program,  statistical  analysis  of  this  information  is  re¬ 
alized  by  making  the  following  two  decisions:  1)  storage  of  the  landscape 
element  features  in  the  computer  memory  for  accumulation  of  a  large  amount 
of  information  ("candidate");  2)  decision-making  on  the  presence  of  anomaly 
and  transmission  of  the  anomaly  information  to  the  ground  facilities. 

The  effectivity  of  the  search  system  operation  will  depend,  in  particu¬ 
lar,  on  how  we  construct  the  selection  algorithm  for  the  "suspected"  landsca¬ 
pe  elements:  on  the  one  hand,  the  storage  of  these  elements  in  the  computer 
memory  decreases  its  oapacity  for  preserving  current  in-line  data,  on  the 
other  hand,  roughening  of  the  criteria  of  selection  of  suspected  elements 
deminishes  the  detection  probability.  The  solution  of  this  contradiction 
lies  in  applying  the  procedure  of  time  exhaustive  search  for  values  of  sepa¬ 
rate  T  components  as  suggested  in  [l]  .  Let  us  consider  the  peculiarities  of 
applying  this  procedure  in  problems  of  global  detection  of  anomlies  on  the 
Earth's  surface  under  monitoring  conditions. 

Let  the  entry  of  the  monitoring  system  have  the  information  going 
through  the  n 1  channels  in  the  form  of  vector  T»(t^,...,tm).  In  this  oaae, 
the  parallel  successive  procedure  for  estimation  of  value  T-T*  based  on  the 
successive  fixation  of  signal  values  of  components  with  the  help  of  m  suc¬ 
cessive  computers  is  used.  The  simultaneous  operation  of  all  m  selection 
devioes  is  achieved  by  matching  the  influxes  of  variants  of  T  component 
values.  ’.Vhile  the  accepted  variants  are  transferred  from  P ^  to 
the  last  has  time  to  process  the  previous  information.  After  Fm,  the 
computer  makes  a  final  estimation  for  T  with  the  speed  V  op/sec.  Taking 
into  account  the  computer  speed  V  of  the  monitoring  system  at  given  period 
of  time  't  for  estimation  of  T  signal  value  we  «  express  the  condition  of 
information  influx  matching  in  terms  of  the  following  formula: 

Vr-tiiA/i-  <*Lr\%NLN%** . . .  •  4,  h**  NL. ..  ¥  {1) 
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where  le  possible  number  of  different  values  of  tit  nt  is  sample  volume; 

<>4  is  probability  of  classifies ting  tA  as  a  signal  one. 

The  presence  of  random  deviations  from  matching  conditions  in  T  in  flux 
requires  the  buffer  memory  for  storage  of  variants  delaying  from  F.  1  in  each 
device. 

Two  types  of  memory  characterized  by  different  forms  of  the  information 
transfer  from  to  F  are  given: 

1)  preserving  the  constant  variant  number; 

2)  preserving  variants  for  given  delay  time. 

In  work  [ 2 ]  were  obtained  the  estimate  of  probability  of  an  error  system 
and  wanted  memory  capacity  with  constant  waiting  time  as  well  as  an  error 
system  and  wanted  delay  memory  capacity  with  constant  number  of  memory  sells. 
The  principle  scheme  considered  above  makes  possible  a  complex  using  the  dif¬ 
ferent  wave  ranges  (optical,  infrared  and  radiorange)  for  the  Earth's  surface 
sensing.  However,  the  problem  concerning  the  realization  of  such  scheme  re¬ 
quires  the  more  comprehensive  consideration. 

The  criteria  of  estimation  of  the  monitoring  system  effectivity  lies  in 
probability  of  solving  of  a  global  problem  of  detection  and  classification  of 
phenomena  on  the  Earth's  surface.  This  probability  depends  on  all  complexes 
of  parameters  of  a  system  and,  in  particular,  on  distribution  of  anomal  sig¬ 
nals  and  noises  (phone). 

One  of  the  possible  models  of  the  monitoring  system  may  be  founded  on  the 
conception  of  "spotness"  of  space.  As  usual,  the  moving  anomaly  reserves  a 
trace  for  itself  whioh  may  change  the  spot  structures  and  sizes.  The  inmov- 
able  anomaly  may  also  cause  the  similar  local  changes  and  be  the  smaller  size 
The  character  of  these  changes  depends  on  the  relation  between  spot  and 
anomaly  sizes.  In  order  to  solve  the  problem  of  detecting  the  anomalies  based 
on  using  "spots",  it  is  necessary  to  have  their  distribution  which,  obviously, 
may  be  obtained  only  by  empirical  way. 

Under  real  conditions,  the  study  of  "spots",  acquisition  of  statistical 
data  on  them  and  the  use  of  this  data  in  a  detector  present  difficulties. 

The  conventional  method  of  detecting  spots  is  the  out-off  selection  method. 
Because  of  this,  the  spot  is  considered  as  a  part  of  space  for  which  the 
parameter  of  environment  for  given  channel  exceeds  or  does  not  exceed  a  cut¬ 
off  value.  The  development  of  algorithm  of  spot  contour  description  and 
creation  of  the  computer  program  for  calculating  their  statistic  parameters 
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wait  for  further  consideration 


The  problems  of  detecting  the  anomalies  are  included  into  the  general 
problem  of  mapping  according  to  data  on  microwave  sensing  the  Earth's  surface 
The  recognition  of  homogeneous  formation  contours  and  the  determination  of 
geophysical  characteristics  may  be  carried  out  by  microwave  radiometry. 
Therefore  the  development  of  new  methods  of  automatic  data  processing  from 
satellite  and  fixing  the  measurements  on  the  surface  are  required. 


After  satellite  K  flights,  we  obt.iln  the  sample  from  K  -  J  processes 
for  two-dimensional  processing.  Here  J  flights  are  excluded  from  studing 
as  no  having  useful  information.  The  knowledge  of  each  flight  course  using 
the  two-dimensional  correlation  analysis  and  interpolation  method  permit  us 
to  make  and  print  a  map  of  radiobrightnass  temperature  which  may  be  used  di¬ 
rectly  when  observing  different  phenomena  on  the  Earth' s  surface.  The  simp- 
liest  algorithm  of  making  such  map  by  a  computer  is  determined  in  the  fol¬ 
lowing  way. 


>*ti)  .a) 
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Wjf  is  reper  object 


Let  7 

be  a  vector  of  registration  data  of  i  flight,  where 

response;  is  flight  time  moment  over  reper  object;  <5 is  angle 

between  course  and  given  constant  direction;  W^1'  are  measurement  results 


After  K  flights,  the  computer  obtains  matrix  //  Wlj  the  lines  of  which 
are  vectors  W™  and  reduces  it  by  removing  ^  lines  as  no  having  the 
response  from  reper  object  or  useful  information. 


It  should  be  noted  that  for  further  operation  jl  ^jl  lines  are  reooried 
so  that  the  V/M  corresponds  to  in  one  direction.  Then  on  placing 

jlWH  lines  in  increasing  order  of  ,  the  brightness  temperature  at  the 

arbitrary  point  M  of  the  region  being  studied  may  be  calculated  ac¬ 

cording  to  the  fonuula 


where  and  ^  are  the  closest  to  . 

There  are  also  other  algorithms  of  interpolation. 

The  development  of  operative  methods  for  determining  the  atmosphere  and 
underlying  surface  geophysical  characteristics  is  considered  to  be  the  impor¬ 
tant  aspect  for  construction  of  the  monitoring  system.  The  automatization  of 
satellite  measurement  processing  involves  the  increase  of  quality  of  mathema¬ 
tical  methods  applied.  The  filtration  of  errors  on  registration  the  telemetric 
signal  and  the  estimation  of  accuracy  of  geophysical  parameter  measurements 
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arc  of  great  Importance*  The  highest  accuracy  of  geophysical  parameter  deter¬ 
mination  may  be  achieved  by  .loint  using  the  information  containing  in  remote 
measurements  and  a  priori  information  which  may  be  obtained  by  other  ways. 

Some  formalization  possibilities  of  remote  sensing  inverse  problem  will 
be  demonstrated  by  way  of  example,  by  processing  microwave  radiometric  measu¬ 
rements  from  the  oceanographic  satellite  Cosmos-1151  [5,6]  •  In  January  1980 
this  satellite  was  launched  into  a  circular  orbit  at  the  altitude  of  650  km 
and  a  slope  angle  of  82,5°  to  equator  plane.  Pour  channel  microwave  radiometer 
mounted  on  the  satellite,  recorded  nadir  radiothermal  radiation  at  the  wave¬ 
lengths  of  0,8;  1,15;  3,2;  8,5  cm.  In  the  course  of  present  experiment  there 
was  developed  the  statistic  method  for  determination  of  atmospheric  water 
vapor  abundance  Q,  liquid  water  content  of  clouds  VV,  ocean  surface  temperatu- 
re  TB  and  wind  speed  V  at  the  oceans  surface  from  microwave  measurements  in 
above  mentioned  channels.  The  method  takes  into  account  statistic  properties 
of  errors  in  absolute  microwave  measurements  and  available  a  priori  informa¬ 
tion  on  the  sought  parameters. 

A  block  diagram  for  automatic  processing  of  satellite  measurements  is 
given  in  Fig.1.  At  the  primary  processing  stage  the  file  of  telemetric  read¬ 
outs  of  the  output  signal  in  operating  channels  of  the  microwave  radiometer 
is  into  a  computer.  The  values  of  the  signal  normalized  to  internal  calibra¬ 
tion  are  determined.  These  normalized  values  are  transferred  to  the  scale  of 
absolute  brightness  temperatures  Tbl>.  For  this  purpose  the  homogeneous  earth 
surface  regions  in  cloudness  conditions  are  used  as  standard  radiators. 

In  parallel  to  Tbr  determination,  the  covariance  matrix  C  of  absolute  measu¬ 
rement  errors  in  various  microwave  channels  is  estimated  by  the  proceedings 
presented  in  [6]  ,  The  obtained  brightness  tempera tures  are  averaged  accord¬ 
ing  to  the  given  space  resolution,  which  partly  makes  up  for  the  difference 
in  antenna  patterns  in  various  channels. 

At  the  secondary  processing  stage  n  geophysical  parameters  along  the 
satellite  orbit  projection  are  determined  from  the  obtained  spectra  of  Tbr. 

The  parameters  are  estimated  by  solving  the  set  of  equations  which  presents 
the  relation  between  geophysical  parameters  and  Tbr  in  m  microwave  radiomet¬ 
ric  channels: 

v  (o)-T^  st  („ 

where  Y(  &J  is  m  -  dimensional  vector  function  of  n  -  dimensional  vector  of 
geophysical  parameters  ;  T  is  m  —  dimensional  vector  of  the  measured  value 
Tbr*  T  is  m  -  dimensional  vector  of  measurement  errors.  In  solution  of  the 
inverse  problem  are  used  additional  limitations  to  determinable  parameters 
taking  into  account  a  priori  information  on  permissible  variation  boundaries 
of  these  parameters.  This  makes  it  possible  to  avoid  senseless  values, 
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involved  with  application  of,  for  example,  the  "quasistical  method"  (negative 
values  of  the  cloud  liquid  water  content  in  Fig, 1,2  £ 7 ]  ). 

As  a  result,  the  optimal  value  ^  is  determined  by  solving  the  ext¬ 
remum  problem 

0* •  arg  min  ^  W  (4) 

ee.F 

where  F  is  the  set  singled  out  in  n  -  space  of  parameters  by  a  priori  limita¬ 
tions,  The  objective  function  ia  constructed  by  the  least-squares 

technique.  In  this  case  the  value  of  the  covariance  error  matrix,  obtained  at 
the  primary  processing  stage,  is  used: 

T(e,  t)  -  rw  -  rjTc"[rw  -r]  (5) 

This  approach  enables  the  information  provided  by  remote  measurements  and  a 
priori  data  on  the  sought  parameters  to  be  optimally  considered  in  the  solu¬ 
tion.  The  least-squares  technique  makes  it  possible  to  estimate  the  accuracy 
of  geophysical  parameter  determination  with  due  regard  for  noise  and  calibra¬ 
tion  measurement  errors.  The  extremum  problem  (4)  is  solved  by  gradient 
methods.  It  does  not  take  much  machine  time  due  to  Bmall  dimension  of  the 
problem.  When  processing  every  following  point  of  measurements,  the  solution 
obtained  for  the  preceding  point  is  used  as  sero-order  approximation.  This 
saves  machine  time  because  the  values  of  geophysical  parameters  in  neigh¬ 
bouring  geographlo  points  stand  close  together. 

The  present  method  has  the  advantage  of  single-valued  relation  between 
the  function  of  <Pf©> 7")  and  the  density  of  posterior  parameter  distribution 
according  to  the  problem  (4).  The  value  of  ¥(0^)  at  its  local  minimum 
point  at  F  may  serve  as  a  quality  criterion  for  the  obtained  solution  in 
automatic  processing.  A  very  large  value  of  eP(^>7~)  indicates  either  a  hard- 
wave  error  or  Inadequacy  of  the  accepted  model  (3).  The  condition 
may  serve  as  a  solution  quality  criterion  [q]  . 

The  meansquare  errors  in  geophysical  parameter  determination  for  the 

2  9 

present  experiment  are  the  following:  Q-0,3  g/cm  ;  W-0,07*0,1  og/m  ; 

V-3+4  m/sec.  The  values  of  Ts  may  be  determined  with  reasonable  accuracy  from 
measurements  in  radiometer  nadir  channels  of  satellite  Cosmos-1151  only  in  the 
absense  of  foaming  caused  by  storm.  Under  these  conditions  the  mean-square 
error  of  Ts  determination  constitites  2-3K. 

In  Fig. 2  is  shown  a  latitudal  distribution  of  Q,\V,  A  V  and  Ts  determined 
from  measurements  with  satellite  Cosmos-1151  over  the  northern  part  of  the 
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Atlantic  Ocean  on  January  25,  1980.  In  conformity  with  available  models,  the 
value  of  d  V  stands  here  for  the  difference  between  V  and  the  critical 
wind  speed  of  V0(V0  &  7  ♦  10  m/sec)  which  induces  foaming.  At  the  back¬ 
ground  of  latitudal  humidity  distribution  one  can  clearly  see  anomalous  re¬ 
gions  of  Q  increase  caused  by  cloudv  systems  of  the  intertropical  convergence 
zone  (2-6  S,L)  and  powerful  atmospheric  front  at  15-50  S,L.  Marked  anticor¬ 
relation  of  V/  and  A  V  values  in  this  front  zone  is  due  to  incomplete  dis¬ 
tinction  between  these  parameter  variations  in  remote  measurements. 

For  comparison,  in  Fig.  3  are  given  the  processing  results  of  the  sance 
measurements  neglecting  a  priori  limitations  to  parameters.  It  can  be  seen 
that  measurement  errors  in  this  case  renult  in  the  negative  values  of  W  and 
A  V  parameters  which  present  difficulties  for  interpretation.  The  negative 
values  of  A  V  obtained  at  2°-4°S,L  and  6°-10°S,L  most  probably  indicate  the 
absence  of  foaming  in  these  regions.  In  Fig. 2  are  given  the  values  of  Ts 
obtained  with  regard  for  AV-0  limitations.  These  values  are  in  good  agre¬ 
ement  with  climatic  data  for  these  regions. 
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Fig. 1 .  Automatic  processing  block  diagram  for 
microwave  radiometric  measurements  from 
Cosmos-1151 . 
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Fig. 2.  Latitudal  distributions  of  Q,,V,  V  and  Ts 

obtained  by  processing  satellite  measurements 
over  northern  part  of  Atlantic  Ocean  on 
January  25,  1980. 


Fig. 3.  Latitudal  distributions  of  Q,w,  V  obtained  by 
processing  the  same  measurements  neglecting  a 
priori  information. 
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